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In SU(3) lattice QCD, we perform the detailed study for the ground-state three-quark (3Q) potential V^q' 
and the 1st excited-state 3Q potential V:jQ ' , i.e., the energies of the ground state and the 1st excited state of the 
gluon field in the presence of the static three quarks. From the accurate calculation for more than 300 different 
patterns of 3Q systems, the static ground-state 3Q potential V^q' is found to be well described by the Coulomb 
plus Y-type linear potential (Y-Ansatz) within 1%-level deviation. As a clear evidence for Y-Ansatz, Y-type flux- 
tube formation is actually observed on the lattice in maximally- Abelian projected QCD. For about 100 patterns 
of 3Q systems, we calculate the 1st excited-state 3Q potential V^q ' , and find a large gluonic-excitation energy 
AEsQ = V'^q' — V^q' of about 1 GeV, which gives a physical reason of the success of the quark model even 
without gluonic excitations. We present also the first study for the penta-quark potential Vsq in lattice QCD, 
and find that Vsq is well described by the sum of the OGE Coulomb plus multi-Y type linear potential. 



1. Introduction 

Quantum chromodynamics (QCD), the SU(3) 
gauge theory, was first proposed by Yoichiro 
Nambu [J in 1966 as a candidate for the funda- 
mental theory of the strong interaction, just af- 
ter the introduction of the "new" quantum num- 
ber, "color" [2]. In spite of its simple form, QCD 
creates thousands of hadrons and leads to vari- 
ous interesting nonperturbative phenomena such 
as color confinement 3 and dynamical chiral- 
symmetry breaking 0]. Even now, it is very diffi- 
cult to deal with QCD due to its strong-coupling 
nature in the infrared region. 

In recent years, the lattice QCD Monte Carlo 
calculation becomes a reliable and useful method 
for the analysis of nonperturbative QCD [S], 
which indicates an important direction in the 
hadron physics. In this paper, using lattice 
QCD, we study the inter-quark potential in de- 
tail |6I7I8I9| . 

In general, the three-body force is regarded as 
a residual interaction in most fields in physics. In 
QCD, however, the three-body force among three 
quarks is a "primary" force reflecting the SU(3) 
gauge symmetry. In fact, the three-quark (3Q) 
potential is directly responsible for the structure 
and properties of baryons, similar to the relevant 



role of the Q-Q potential for meson properties, 
and both the Q-Q potential and the 3Q poten- 
tial are equally important fundamental quanti- 
ties in QCD. Furthermore, the 3Q potential is 
the key quantity to clarify the quark confinement 
in baryons. However, in contrast to the Q-Q po- 
tential [Sj , there was almost no lattice QCD study 
for the 3Q potential before our study in 1999 
in spite of its importance in the hadron physics. 

2. The Ground-State 3Q Potential in QCD 

The Q-Q potential is known to be well de- 
scribed with the inter-quark distance r as |5I6I7 | 

^qqM = -^ + W + ^qq- (1) 

As for the 3Q potential form, we note two theoret- 
ical arguments at short and long distance limits. 

1. At the short distance, perturbative QCD is ap- 
plicable, and therefore 3Q potential is expressed 
as the sum of the two-body Coulomb potential 
originating from the onc-giuon-cxchange process. 

2. At the long distance, the strong-coupling ex- 
pansion of QCD is plausible, and it leads to the 
flux-tube picture ^O!- For the 3Q system, there 
appears a junction which connects the three flux- 
tubes from the three quarks, and Y-type flux-tube 
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formation is expected |6I7I1U| . 

Then, we theoretically conjecture the func- 
tional form of the 3Q potential as the Coulomb 
plus Y-type linear potential, i.e., Y-Ansatz, 

^3q'' = ~^3Q Ij.. + '^SQLmin + C3Q, (2) 

i<j ' ^ ' 

where Lmin is the minimal value of the total flux- 
tube length. Of course, it is nontrivial that these 
simple arguments on UV and IR limits of QCD 
hold for the intermediate region. Then, we study 
the 3Q potential in lattice QCD. Note that the 
lattice QCD calculation is completely indepen- 
dent of any Ansatz for the potential form. 

2.1. The Three-Quark Wilson Loop 

Similar to the Q-Q potential calculated with 
the Wilson loop, the 3Q potential can be calcu- 
lated with the 3Q Wilson loop |6I7I8| defined as 

W'sQ = ^e.tcea'b'c'Ur'ufuf (3) 

with Uk = Pexp{ig J^^ dXf,A^'{x)} in Fig.l. The 
3Q Wilson loop physically means that a color- 
singlet gauge-invariant 3Q state is created at 
t — and is annihilated at t — T with the three 
quarks spatially fixed for < t < T. 




Figure 1. The 3Q Wilson loop W3Q. 

The vacuum expectation value of the 3Q Wil- 
son loop is expressed as 

00 

n=0 

where Vn denotes the nth energy of the gauge- 
field configuration in the presence of the spatially- 
fixed three quarks |6I7I8| . 

While Vn depends only on the 3Q location, C„ 
depends on the operator choice for the 3Q state. 



For the accurate calculation, we need to reduce 
the excited-state component, i.e., C„(n > 1), in 
the 3Q state prepared at t = and T. 

2.2. The Smearing Method 

The smearing is a useful method to construct 
the quasi-ground-state operator in lattice QCD 
in a gauge-invariant manner, and is actually suc- 
cessful for the ground-state Q-Q potential [Jj. 
(Note that the smearing is just a method to con- 
struct the operator, and hence it never changes 
the gauge configuration unlike the cooling.) 

The smeared link- variable C/^(s) includes a spa- 
tial extension, and the smeared "line" expressed 
with Ufj^{s) corresponds to a Gaussian-distributed 
"flux-tube" in terms of the original link- variable 
[7]. Therefore, the properly smeared line is ex- 
pected to resemble the ground-state flux-tube. 

Then, through the selection of the properly 
smeared 3Q Wilson loop (W3q[C^^(s)]), we can 
construct the ground-state-dominant 3Q operator 
for the accurate measurement of the ground-state 
3Q potential 151715] . 

2.3. Lattice QCD results for 3Q Potential 

For more than 300 different patterns of 
spatially-fixed 3Q systems, we perform the thor- 
ough calculation of the ground-state potential 
V^Q in SU(3) lattice QCD with the standard 
plaquette action with 12'^ x 24 at /3 = 5.7 and 
with 16^ X 32 at /3=5.8 and 6.0 at the quenched 
level. For the accurate measurement, we use the 
smearing method and construct the ground-state- 
dominant 3Q operator |()l7l8j . 

To conclude, we find that the static ground- 
state 3Q potential V^q' is well described by the 
Coulomb plus Y-type linear potential (Y-Ansatz) 
within 1%-level deviation |HI7| . 

2.4. Other Studies on the 3Q Potential 

To clarify the current status of the 3Q poten- 
tial, we introduce recent works of other groups. 

de Forcrand's group, who once supported A- 
Ansatz in lattice QCD 11 , seems to change their 
opinion from A- Ansatz to Y-Ansatz |12| . 

Kuzmenko and Simonov also showed that the 
Delta-shape is impossible from gauge-invariance 
point of view, and the Y-shaped configuration is 
the only possible for the three-quark system ^3] . 
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One of the theoretical basis of A-Ansatz was 
CornwaU's conjecture based on the vortex vac- 
uum model |14| . Recently, Cornwall re-examined 
his previous work and found that the correct an- 
swer is Y-Ansatz instead of A-Ansatz |15j . 

In this way, Y-Ansatz for the static 3Q poten- 
tial seems almost settled both in lattice QCD and 
in analytic framework. 

2.5. Y-type Flux- Tube Formation 

Recently, as a clear evidence for Y-Ansatz, Y- 
type flux-tube formation is actually observed in 
maximally- Abelian (MA) projected lattice QCD 
from the measurement of the action density in the 
spatially-fixed 3Q system I6"T7" . (See Fig. 2.) 




Figure 2. The lattice QCD result for Y-type flux- 
tube formation in the spatially-fixed 3Q system 
in MA projected QCD. The distance between the 
junction and each quark is about 0.5 fm. 

3. Penta-Quark Potential in Lattice QCD 

Motivated by the recent discovery of the penta- 
quark baryon O+(1540), we perform the first 
study of the static penta-quark (5Q) potential 
VsQ in SU(3) lattice QCD with /3=6.0 and le^x 
32 at the quenched level. We investigate the QQ- 
Q-QQ configuration as shown in Fig. 3. With the 
smearing method j6l7l8j to enhance the ground- 
state component, we accurately calculate the 5Q 
potential V5Q from the 5Q Wilson loop (Wsq) as 
shown in Fig. 4 in a gauge-invariant manner. 

We find that the 5Q potential V5Q is well de- 
scribed by the sum of one-gluon-exchange (OGE) 
Coulomb term and multi-Y type linear term, 
which wc call the "OGE plus multi-Y Ansatz" , 
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Figure 3. (a) The QQ-Q-QQ type configuration 
for the penta-quark system. The two QQ clus- 
ters belong to the 3* representation of the color 
SU(3). (b) A planar configuration of the penta- 
quark system. We take di — d2 = ds = di = d. 
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Figure 4. The penta-quark (5Q) Wilson loop 
W5Q for the 5Q potential V5Q. The contours 
M, Li, Ri{i — 3,4) are line-like and Li,Ri{i = 
1, 2) are staple-like. The 5Q gauge- invariant state 
is generated at i = and is annihilated at t = T. 
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with r, 



r^l and Lmin 



r45 



being the minimal 
length of the flux-tube linking five quarks. Note 
that there appear three kinds of Coulomb coef- 
ficients (^5Q, ^^5Q, j^5q) in the penta-quark 
system, while only one Coulomb coefficient, Aqq 
or A3Q, appears in the QQ or the 3Q system. 
(In Eq.(|Hll, A5Q corresponds to A^q or ^^qq in 
terms of the OGE result.) 

Figure 5 shows the lattice QCD results for the 
5Q potential V5Q. The symbols denote the lattice 
data, and the curves denote the theoretical form 
of the OGE plus multi-Y Ansatz with (A5Q,cr5Q) 
fixed to be (A3q,CT3q) in the 3Q potential V3Q in 
Ref.[7]. (Note that there is no adjustable parame- 
ter for the theoretical curves besides an irrelevant 
constant C5Q, since A5Q and ctsq are fixed.) One 
finds a good agreement between the lattice data 
and the OGE plus multi-Y Ansatz. 

We note that the multi-quark system including 
four or more quarks can take a three-dimensional 
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Figure 5. Lattice QCD results of the penta- 
quark potential V5Q for the planar 5Q configura- 
tion with hi = h2 = h in Fig. 3 in the lattice unit. 
Each 5Q system is labeled by d and h. The sym- 
bols denote the lattice data, and the curves the 
theoretical form of the OGE plus multi-Y Ansatz. 

shape, while the QQ and the 3Q systems can take 
only planar configuration. Then, we investigate 
also the twisted 5Q configuration, and find that 
the planar and the twisted 5Q configurations are 
almost degenerate. Then, no special configura- 
tion is favored in the 5Q system in terms of the 
energy, and general 5Q systems tend to take a 
three-dimensional configuration. 

From the comparison with the QQ and the 3Q 
potentials 6 7 8, the universality of the string 
tension and the OGE result are found among QQ, 
3Q and 5Q systems as 

f^QQ - ^3Q - -Aqq ~ ^3Q ~ AsQ. (6) 

This result supports the flux-tube picture on the 
confinement mechanism even for the multi-quark 
system. 

4. Gluonic Excitations in the 3Q System 

In 1969, Y. Nambu first pointed out the string 
picture for hadrons jl8| . Since then, the string 
picture has been one of the most important pic- 
tures for hadrons and has provided many interest- 
ing ideas in the wide region of the particle physics. 

For instance, the hadronic string creates in- 
finite number of hadron resonances as the vi- 
brational modes, and these excitations lead to 
the Hagedorn "ultimate" temperature 19 , which 
gives an interesting theoretical picture for the 



QCD phase transition. 

For the real hadrons, of course, the hadronic 
string is to have a spatial extension like the flux- 
tube, as the result of one-dimensional squeez- 
ing of the color-electric flux in accordance with 
color confinement Therefore, the vibrational 
modes of the hadronic flux-tube should be much 
more complicated, and the analysis of the excita- 
tion modes is important to clarify the underlying 
picture for real hadrons. 

In the language of QCD, such non-quark-origin 
excitation is called as the "gluonic excitation", 
and is physically interpreted as the excitation of 
the gluon-field configuration in the presence of 
the quark-antiquark pair or the three quarks. 

In the hadron physics, the gluonic excitation 
is one of the interesting phenomena beyond the 
quark model, and relates to the hybrid hadrons 
such as qqG and qqqG. In particular, the hy- 
brid meson includes the exotic hadrons with 
J^^ = 0--,0+-,l-+,2+-,---, which cannot be 
constructed within the simple quark model. 

In this section, we study the excited-state 3Q 
potential and the gluonic excitation using lattice 
QCD |8j, to get deeper insight on these exci- 
tations beyond the hypothetical models such as 
the string and the flux-tube models. Here, the 
excited-state 3Q potential is the energy of the 
excited state of the gluon-field configuration in 
the presence of the static three quarks, and the 
gluonic-excitation energy is expressed as the en- 
ergy difference between the ground-state 3Q po- 
tential and the excited-state 3Q potential. 

4.1. General Formalism 

We present the formalism to extract the 
excited-state potential 0. For the simple nota- 
tion, the ground state is regarded as the "0th ex- 
cited state". For the physical eigenstates of the 
QCD Hamiltonian H for the spatially-fixed 3Q 
system, we denote the nth excited state by \n) 
[n = 0, 1, 2, • ■ •). Since the three quarks are spa- 
tially fixed in this case, the eigenvalue of H is 
expressed by a static potential as H\n) = Vn\n), 
where Vn denotes the nth excited-state 3Q po- 
tential. Note that both Vn and jn) are universal 
physical quantities relating to the QCD Hamil- 
tonian H. In fact, Vn depends only on the 3Q 
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location, and \n) satisfies the orthogonal condi- 
tion as {mln) = (5„m- 

Suppose that (fc = 0, 1, 2, • ■ •) are arbitrary 
given independent spatially-fixed 3Q states. In 
general, each 3Q state can be expressed by a 
linear combination of 3Q physical eigenstates \n), 



\^k)^4\o)+4\i)+4m + ' 



(7) 



Here, the coefficients cjij depend on the selection of 
|$fc), and hence they are not universal quantities. 

The Euclidean-time evolution of the 3Q state 
$fc(t)) is expressed with the operator e~^*, 
which corresponds to the transfer matrix in lat- 
tice QCD. The overlap ($j(T)|$fe(0)) is given by 
the 3Q Wilson loop with the initial state |$fc) at 
t = and the final state at t = T, and is 
expressed in the Euclidean Heisenberg picture as 

W^' ^ ($,|iy3Q(r)|<i>fc) = (<i>,|e-^^|<i>,) 



ni—O n—Q n—Q 

Using the matrix C satisfying C"'' — cjj and the 
diagonal matrix as A™" = e^^"'^S™^, we 
rewrite the above relation as 



(9) 



Note here that C is not a unitary matrix, and 
hence this relation does not mean the simple di- 
agonalization by the unitary transformation. 

Since we are interested in the 3Q potential Vn 
in At rather than the non-universal matrix C, we 
single out Vn from the 3Q Wilson loop Wt as 



VF-iR^T+i = C"Miag(e" 



1 1 



(10) 



which is a similarity transformation. Then, e~^" 
can be obtained as the eigenvalues of the matrix 
Wt^Wt+1, i.e., solutions of the secular equation. 



det{W^^WT+i - tl} = J|(e 



t)-0. (11) 



Thus, the 3Q potential Vn can be obtained from 
the matrix Wt^Wt+i- 

In the practical calculation, we prepare TV inde- 
pendent sample states \^k) {k — 0,1, ■ ■ ■ , N — 1). 
By choosing appropriate states \^k) so as not 
to include highly excited-state components, the 



physical states \n) can be truncated as < n < 
iV — 1. Then, Wt, C and are reduced into 
N X N matrices, and the secular equation Hll|l 
becomes the TVth order equation. 

4.2. Lattice QCD for Gluonic Excitations 

For about 100 different patterns of spatially- 
fixed 3Q systems, we calculate the excited-state 
potential using SU(3) lattice QCD with 16^ x 32 
at (3=5.8 and 6.0 at the quenched level 8 . In 
Fig. 6, we show the 1st excited-state 3Q potential 



Vg^Q ■ and the ground-state potential V^q ' . 




5 10 
L„i„ [lattice unit] 



e.s. 

3Q 



Figure 6. The 1st excited-state 3Q potential V^ 
and the ground-state 3Q potential V^q' . The lat- 
tice results at /3 = 5.8 and (3 = 6.0 well coincide 
besides an irrelevant overall constant. The glu- 
onic excitation energy AE3Q = V^q- — V^q' is 
found to be about IGeV in the hadronic scale. 



The energy gap between V^q' and V^q - phys- 
ically means the excitation energy of the gluon- 
field configuration in the presence of the spatially- 
fixed three quarks, and the gluonic excitation en- 
ergy AE3Q = V^Q - — V^i^' is found to be about 
IGeV |8I17| in the hadronic scale as Lmin ~ 1 fm. 

Note that the gluonic excitation energy of 
about IGeV is rather large in comparison with 
the excitation energies of the quark origin. The 
present result predicts that the lowest hybrid 
baryon, which is described as qqqG in the va- 
lence picture, has a large mass of about 2 GeV 
|8ll7j . (The present result seems to suggest the 
constituent gluon mass of about IGeV.) 

5. Behind the Success of the Quark Model 

Finally, we consider the connection between 
QCD and the quark model in terms of the gluonic 
excitation |8I17| . While QCD is described with 
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quarks and gluons, the simple quark model suc- 
cessfully describes low-lying hadrons even with- 
out explicit gluonic modes. In fact, the gluonic 
excitation seems invisible in the low-lying hadron 
spectra, which is rather mysterious. 

On this point, we find the gluonic-excitation 
energy to be about IGeV or more, which is rather 
large compared with the excitation energies of the 
quark origin, and therefore the effect of gluonic 
excitations is negligible and quark degrees of free- 
dom plays the dominant role in low-lying hadrons 
with the excitation energy below IGeV. 

Thus, the large gluonic-excitation energy of 
about IGeV gives the physical reason for the in- 
visible gluonic excitation in low- lying hadrons, 
which would play the key role for the success of 
the quark model without gluonic modes )8ll7j . 

In Fig. 7, we present a possible scenario 
from QCD to the massive quark model in 
terms of color confinement and dynamical chiral- 
symmetry breaking (DCSB) |17|. 
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Quantum Chromodynamics 
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Figure 7. A possible scenario from QCD to the 
quark model in terms of color confinement and 
DCSB. DCSB provides a large constituent quark 
mass of about 300McV, which enables the non- 
relativistic treatment for quark dynamics. Color 
confinement leads to the color flux-tube forma- 
tion among quarks. In the flux-tube picture, the 
gluonic excitation is described as the flux-tube 
vibration, and then its energy is expected to be 
large in the hadronic scale. The large gluonic- 
excitation energy of about IGeV leads to the ab- 
sence of the gluonic mode in low-lying hadrons, 
which would play the key role to the success of the 
quark model without gluonic excitation modes. 



